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Abstract 

The B meson lifetime ratios are calculated to the order of in the 

heavy quark expansion. The predictions of those ratios are dependent on 
four unknown hadronic parameters -Bi, B2, £1 and £2-, where Bi and B2 
parametrize the matrix elements of color singlet-singlet four-quark opera- 
tors and £1 and £2 the matrix elements of color octet-octet operators. We 
derive the renormalization-group improved QCD sum rules for these param- 
eters within the framework of heavy quark effective theory. The results are 
5i(mb) = 0.96 ± 0.04, B2{rrn,) = 0.95 ± 0.02, ei(mfc) = -0.14 ± 0.01, and 
£2{mb) = —0.08 lb 0.01 to zeroth order in 1/mf,. The resultant B meson life- 
time ratios are T{B-)/T{Bd) = 1.11 ± 0.02 and T{Bs)/T{Bd) ^ 1 in SU(3) 
symmetry limit. 
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I. INTRODUCTION 



A QCD-based formulation for treatment of inclusive heavy hadron decays has been de- 
veloped in past years [1-3]. According to the optical theorem, the inclusive decay rates 
are related to the imaginary part of certain forward scattering amplitudes along the physi- 
cal cut. Since the cut is dominated by physical intermediate hadron states like resonances 
which are nonperturbative in nature, a priori the operator product expansion (OPE) or 
heavy quark expansion cannot be carried out on the physical cut. Nevertheless, for inclusive 
semileptonic decays, OPE can be employed for some smeared or averaged physical quanti- 
ties. For example, by integrating out the neutrino energy, one can apply the OPE to the 
double differential cross section d?T /{dcf'dEf) by deforming the contour of integration into 
the unphysical region far away from the physical cut [4]. Therefore, global quark-hadron 
duality, namely the matching between the hadronic and OPE-based expressions for decay 
widths and smeared spectra in semileptonic B or bottom baryon decays, follows from the 
OPE and is justified except for a small portion of the contour near the physical cut which 
is of order Aqcb/ttiq- Unfortunately, there is no analogous variable to be integrated out in 
inclusive nonleptonic decays, allowing an analytic continuation into the complex plane. As 
a result, one has to invoke the assumption of local quark-hadron duality in order to appply 
the OPE in the physical region [5]. It is obvious that local duality is theoretically less firm 
and secure than global duality. In order to test the validity of local quark-hadron duality, it 
is thus very important to have a reliable estimate of the heavy hadron lifetimes within the 
OPE framework and compare them with experiment. 

In the heavy quark limit, all bottom hadrons have the same lifetimes, a well-known 
result in the parton picture. With the advent of heavy quark effective theory and the OPE 
approach for the analysis of inclusive weak decays, it is realized that the first nonperturbative 
correction to bottom hadron lifetimes starts at order 1/ml and it is model independent (for 
a review, see [6]). However, the 1/ml corrections are small and essentially canceled out 
in the lifetime ratios. The nonspectator effects such as VF-exchange and Pauh interference 
due to four-quark interactions are of order l/niQ, but their contributions can be potentially 
significant due to a phase-space enhancement by a factor of IGtt^. As a result, the lifetime 
differences of heavy hadrons come mainly from the above-mentioned nonspectator effects. 



2 



The world average values for the hfetime ratios of bottom hadrons are [7] : 

t{B- 



r(A.) 



o( = 1.07 ±0.04, 
0.95 ±0.05, 



0.78 ±0.06. (1.1) 



r{Bl) 

Since the model-independent prediction of T{Ah)/T{Bd) to order 1/m^ is very close to unity 
[see Eq. (2.7) below], the conflict between theory and experiment for this lifetime ratio is quite 
striking and has received a lot of attention [8-15]. One possible reason for the discrepancy 
is that local quark-hadron duality may not work in the study of nonleptonic inclusive decay 
widths. Another possibility is that some hadronic matrix elements of four-quark operators 
are probably larger than what naively expected so that the nonspectator effects of order 
IGn^/ml may be large enough to explain the observed lifetime difference between the A;, 
and Bd- Therefore, as stressed by Neubert and Sachrajda [11], one cannot conclude that 
local duality truly fails before a reliable field-theoretical calculation of the four-quark matrix 
elements is obtained. Contrary to the 1/m^ corrections, the estimate of the nonspectator 
effects is, unfortunately, quite model dependent. 

Conventionally, the hadronic matrix elements of four-quark operators are evaluated using 
the factorization approximation for mesons and the quark model for baryons. However, as 
we shall see, nonfactorizable effects absent in the factorization hypothesis can affect the 
B meson lifetime ratios significantly. In order to have a reliable estimate of the hadronic 
parameters Bi, B2, e\ and £2 in the meson sector, to be introduced below, we will apply the 
QCD sum rule to calculate these unknown parameters. After a brief review on the status 
of the OPE approach for the B hadron lifetime ratios in Sec. 11, we proceed to derive in 
Sec. Ill the renormalization-group improved QCD sum rules for the parameters B^ and Si 
and present a detailed analysis. Sec. IV gives discussions and conclusions. 



II. A BRIEF OVERVIEW 

Within the heavy quark expansion framework, we will focus in this paper on the study 
of the four-quark matrix elements of the B meson to understand the problem with B meson 
lifetime ratios. Before proceeding, let us briefly review the content of the theory. Applying 
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the optical theorem, the inclusive decay width of the bottom hadron Hi, containing a b quark 
can be expressed in the form 

r(Hb ^X)^^Im [ d^x {H,\ T{ i£eff(^), >Ceff(0) } \H,) , (2.1) 

where C^q is the relevant effective weak Lagrangian that contributes to the particular final 
state X. When the energy release in a 6 quark decay is sufficiently large, it is possible to 
express the nonlocal operator product in Eq. (2.1) as a series of local operators in powers 
of l/m^ by using the OPE technique. In the OPE series, the only locally gauge invariant 
operator with dimension four, hUph, can be reduced to mi^h by using the equation of motion. 
Therefore, the first nonperturbative correction to the inclusive B hadron decay width starts 
at order 1/m^.^ As a result, the inclusive decay width of a hadron Hi, can be expressed 
as [1,2] 



^ - 1927r3 "3 Wbb\H,) + c, 



^^^x(.)(^y^|H)+o(lK)L (2.2) 



where Vckm denotes some combination of the Cabibbo-Kobayashi-Maskawa parameters and 
cf reflect short-distance dynamics and phase-space corrections. The matrix elements in 
Eq. (2.2) can be systematically expanded in powers of 1/mi, in heavy quark effective theory 
(HQET) [17], in which the 6-quark field is represented by a four- velocity-dependent field 
denoted by h!^^ (x) . 

In Eq. (2.2) cf are functions of Ci and C2, the Wilson coefficients in the effective Hamil- 
tonian 



AB=1 _ 
eff — 



V2 

+ K6K;(ci(/^)0^(/x)+C2(/x)02^(/x)) + ---] +h.c., (2.3) 



where q — d,s, and 



^It is emphasized in [16] that the cancellation of the l/mg corrections to the inclusive decay width 
occurs when it is expressed in terms of the running short-distance quark mass, e.g. the MS mass, 
rather than the pole quark mass. 



or = C7^(l - 75)& - l^)u, = 57^(1 - Ts)^ ^^(1 - Ts)^^ ■ (2.4) 

The scale and scheme dependence of the Wilson coefficients ci,2(a*) are canceled out by the 
corresponding dependence in the matrix element of the four-quark operators Oi,2- That is, 
the four-quark operators in the effective theory have to be renormalized at the same scale 
H and evaluated using the same renormalization scheme as that for the Wilson coefficients. 
Schematically, we can write (Hes) — c{ii){0{ii)) — c{ii)g{ii){0)tree — c^^{0)tree, where the 
effective Wilson coefficients cf^ are renormalization scale and scheme independent. Then the 
factorization approximation or the quark model is applied to evaluate the hadronic matrix 
elements of the operator O at tree level. The explicit expression for g{fi), the perturbative 
corrections to the four-quark operators renormalized at the scale /j,, has been calculated in 
the hterature [20,21]. To the next-to- leading order (NLO) precision, we have [22] ^ 

cf = 1.149 , cf = -0.325 . (2.5) 

Replacing a by cf and using mt, = (4.85 ± 0.25) GeV, {rric/mbf =0.089, \Vcb\ = 0.039, 
Gb^,b, = 0.366 GcV^ Gb, = 0.381 GeV' [23], Ga, = 0, Kb^,b, ~ i^A, ~ 0.4 GeV^, Kb, = 
1.02 X KB^,Bd [24] together with the nonleptonic inclusive results to the next-to- leading order 
[18], we find numerically 

rsL(S ^ ei^X) = (4.18^^:20) X 10-1^ GeV, 
rsL(A6 ^ ePX) = {4.32t\f,) X 10"^^ GeV, 

T{B) = TnLiB) + 2.24rgL(5 ^ ePX) = (3.611^:°^) x 10"^^ GeV, 

r{Ab) = rNL(A6) + 2.24 rsL(A6 ^ ePX) = (3.651^:°^) X lO'^^ Q^y _ (3.6) 



^The effective Wilson coefficients given in (2.5) are derived from Cj(/u) at /U = to the NLO [12]. 
Nevertheless, it is not difficult to explicitly check the scale and scheme independence of cf. For 
example, the authors of [21] obtain cf = 1.160 and cf = —0.334 at /x = 2.5 GeV. Therefore, cf 
are very insensitive to the chosen fi scale, as it should be. It is known that the Wilson coefficient 
C2 at the NLO: C2 = —0.185 in the naive dimension regularization scheme and C2 = —0.228 in 
the 't Hooft-Veltman scheme [20], deviates substantially from the leading-order value C2 = —0.308 
at = mft(m(,). However, the resultant cf is scheme independent and its value is close to the 
leading-order one. 
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It follows that the lifetime ratios of the Hi, hadrons are 

■l + 0{l/ml), 



^^^'^ - 1.0005 + 0(l/mg), 



= 0.99 + C>(l/m^). (2.7) 



"^(^b) _ n no , r^/^ /™3n 



r{Bci) 

Note that t{Bs) here refers to the average lifetime of the two CP eigenstates of the Bs 
meson. It is evident that the corrections are too small to explain the shorter lifetime 

of the Aft relative to that of the B^. To the order of the nonspectator effects due to 

Pauli interference and VF-exchange parametrized in terms of the hadronic parameters [11]: 
Bi, B2, El, 62, B, and r (see below), may contribute significantly to lifetime ratios due to a 
phase-space enhancement by a factor of IGvr^. The four-quark operators relevant to inclusive 
nonleptonic B decays are 

Ol_p = Ir ql Ql bn , 

n-A = bL7^t''qLqLYt''bL, 

Tl^p^bRt^LgLebR, (2.8) 

where qr^l = ^^^Q and = A"/2 with A'' being the Gell-Mann matrices. For the matrix 
elements of these four-quark operators between B hadron states, we follow [11] to adopt the 
following definitions: 

^ {B,\0^y_^\B,)^^-i^B^, 



2mB, 



2mB, 



{Bq\T^_^\Bq)= -El , 



2mB, 



-{Bq\T^_p\Bq)= -82 , 



2mB/ 8 

■(A6|0^_^|Ab)= -T^r, 



2mA/ 48 
^(A.|T^,|A.). -liB + ^)^(A.|0^_.|A.) , (2.9) 



where fs^ is the Bq meson decay constant defined by 

(0|g7^756|S,(p))=i/B^p^. (2.10) 

Under the factorization approximation, Bi — 1 and Ei — 0, and under the valence quark 
approximation B — 1 [11]. 

To the order of 1/ml, we find that the S-hadron hfetime ratios are given by 

^ = 1 + (y^^)'(0.043S, + 0.0006S. - 0.61e, + 0.17s.) , 

= 1 + (t^F^T^) Vl-7 X 10-5 Bi + 1.9 X 10-5 B2 - 0.0044£i + 0.0050 £2) , 

T{B^) ^185 MeV^ 

-0.15£i + 0.17£2 - (0.014 + 0.0195)r] . (2.11) 

The above results are similar to that given in [11]. We see that the coefficients of the 
color singlet-singlet operators are one to two orders of magnitude smaller than those of the 
color octet-octet operators. This implies that even a small deviation from the factorization 
approximation £j = can have a sizable impact on the lifetime ratios. It was argued in 
[11] that the unknown nonfactorizable contributions render it impossible to make reliable 
estimates on the magnitude of the lifetime ratios and even the sign of corrections. That is, 
the theoretical prediction for T{B^)/T{Bci) is not necessarily larger than unity. In the next 
section we will apply the QCD sum rule method to estimate the aforementioned hadronic 
parameters, especially Si. 



III. THE QCD SUM RULE CALCULATION 

In this section we will employ the method of QCD sum rules within the framework of 
HQET. Since the h quark is treated as a static quark with an infinite quark mass in HQET 
and since HQET is a low-energy effective theory, it is natural to regard the matrix elements of 
(2.9) as defined at the scale nib and to evaluate the corresponding hadronic matrix elements 
in HQET at a scale ji < mi,. Indeed, it has been argued that the estimate of hadronic 
matrix elements of four-quark operators using the factorization hypothesis for mesons and 
the quark model for baryons becomes more reliable if the operators are renormalized at a 
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typical hadronic scale [26] . In the sum rule approach, the correlation function (or the so-called 
Green function), within the QCD framework, can be expanded as a series of operators On(Af) 
multiplied by the Wilson coefficients Cn{—'2uj/ix,gs)/oj'^., where u is an external momentum 
flowing in (or out) the correlation function and /x is the factorization scale that separates 
the long-distance part 0„ from the short-distance one C„. The quality of the convergence 
of the OPE series is controlled by the value of the external momentum u. The factorization 
scale jji cannot be chosen too small, otherwise the strong coupling constant would be so 
large that Wilson coefficients cannot be perturbatively calculated. Four-quark operators are 
sometimes renormalized at a typical scale //^ 0.67 GeV, corresponding to the coupling 
constant a^^{iJLh) ~ 0{1). However, such a scale is not quite suitable for the sum rule 
calculation. Instead we choose ji — 1 GeV as the lowest possible factorization scale in 
the ensuing study. After summing over the logarithmic dependence a^XvT' {—2uj / ix) by the 
renormalization- group method, one obtains the nonperturbative quantity X[ij) which can 
be extracted from the correlation function in the following form 

^ C„(l,g,(-2a;)) / «,(//) x7n-E,7. 

where 7^ are the anomalous dimensions of 0„ and 7-, the anomalous dimensions of the cur- 
rents appearing in the correlation function. As noted above, the four-quark operators in 
Eq. (2.9) are defined at the scale // = ruf,. In HQET where the b quark is treated as a static 
quark, we can use the renormalization group equation to express them in terms of the oper- 
ators renormalized at a scale Aqcd <C /x ^ m^. These operators have the hybrid anomalous 
dimensions [27-29] and their rcnormalization-group evolution is determined by the anoma- 
lous dimensions in HQET. The operators and and similarly 0|_p and Tg_p, 
mix under renormalization. In the leading logarithmic approximation, the renormalization- 
group equation of the operator pair (O, T) governed by the hybrid anomalous dimension 
matrix reads ^ 



^One of the off-diagonal anomalous dimension matrix elements in Eq. (3.2) has a sign opposite to 
that obtained in [11], but the final result in Eq. (3.6) is in full agreement with the results derived 
there. 
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where t = | ln((5^//^^), Cp = {N^ — l)/2Nc, and effects of penguin operators induced from 
evolution have been neglected. 

The solution to the evolution equation Eq. (3.2) has the form 



f1 




f 8 2 \ 
9 3 






[t) 


Q 


\ 27 9 / 


\ 


1/ 



D 



n 1 



(3.3) 



where 



D, 



\D2 




(3.4) 



(3.5) 



and I3q 
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■nj is the leading-order expression of the /3-function with Uf being the 



number of light quark flavors. The subscript ^ in Eq. (3.4) and in what follows denotes the 
renormalization point of the operators. Given the evolution equation (3.3) for the four-quark 
operators, we see that the hadronic parameters Bi and Si normalized at the scale mi, are 
related to that aX ji — l GeV by 



( u\ / 8 2 \ / ^9/(2/3o) W Bi-\ei^ 

l)\\Bi + e,j 




> _A 8 . 

\ 27 9 / 



(3.6) 



H=l GeV 



and hence 



Bi{mk) ~ IMBiin) - 0.4l£i(/x) , 
£i(m6) ~ -0.090Si(At) + 1.07£,(At) , 

with a — 1 GeV, where uses have been made of asimz) — 0.118, 
4.85 GeV, = 1.45 GeV and 

ln(lnf^ 



MS 



(3.7) 

333 MeV, = 



a.(Q) 



47r 
/?oln^ 



1 - 



PI Inf 
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to the next-to-leading order with (3i — 51 — ^Uf. The above results (3.7) indicate that 
renormalization effects are quite significant. 

It is easily seen from Eqs. (3.3) and (3.4) that the normalized operator Di (or D2) is 
simply multiplied by Lq^'^^°^ (or 1) when it evolves from a renormalization point to another 
point Q. In what follows, we will apply this property to derive the renormalization-group 
improved QCD sum rules for Dj at the typical scale fj, — 1 GeV. We define the new four-quark 
matrix elements as follows 



1 



B,\Df{i,)\B,)^ 



(3.8) 



where the superscript (i) denotes {V — A) four-quark operators for i = 1 and (5" — P) 
operators for i = 2, and df^ satisfy 



\ "2 



Si 



(3.9) 



Since the terms linear in four-quark matrix elements are already of order we only 

need the relation between the full QCD field h{x) and the HQET field h^^\x) to the zeroth 

order in l/rrif,: h{x) = e~'^"^''^'^ {h^^x) + 0{l/mi,)}. Therefore, in analogue to Eq. (2.8), we 
define the relevant four-quark operators in HQET as 





= h^vht^QL QLl^h^vl , 




^S-P 


— n-vR Ql rivR ' 




rpv 

^V-A 






rpv 




(3.10) 



The corresponding hadronic matrix elements of these four-quark operators are parametrized 

by 

^(^(^)|0^>-Al^(^))^^^^i?r(/^), 
\{B{v)\Ol_AB{v))^^^^Bl{i,\ 

\{B{v)Tv-a\B{v)) ^ ^^elif,), 
l{B{v)\T^s-p\B{v))^^^el{fi), 



(3.11) 
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where the heavy-flavor-independent decay constant F deflned in the heavy quark hmit is 

given by 

{0\qri5h'^^\B{v)) = iF{fx)v'' . (3.12) 

The decay constant -F(/i) depends on the scale // at which the effective current operator is 
renormalized and it is related to the scale-independent decay constant /b of the B meson by 

F{mk) = fB V^. (3.13) 

Notice that F in Eq. (3.11) is chosen to be normalized at the scale m^. 

To complete the aim of obtaining the matrix elements of four-quark operators, we apply 
the method of QCD sum rules [30] . We consider the following three-point correlation function 

<f (a;,a;0 = [ dx dy e^^^-^-'^'^-y {0\T{[q{x)T^h^^\x)] Df\o) [q{y)T ph^^\y)n\0) , (3.14) 



of the operator of^ deflned in Eq. (3.4), where Fq = 7^75. However, this current interpolates 
not only the heavy mesons with quantum number = 0~ but also that with quantum 
number = 1+. Therefore, we need to decompose into Fq, = F^^ — VaT^^, with 
F^^ = (7 + v)a'y5 for — 1"*" and F^*^ = 75 for — 0~. As a consequence, H^^ is recast 
to 



D 



^ {-9af3 + VaVp)Uj^,^i) +VaVf3llj^,(^i), (3.15) 

j 3 

where 



'-9ai3 + VaVj3)Ti'^(i) = dx dy 6 



lujv-x—iLO v-y 



X {0\T{[qix)ri^h^:\x)]Df{0) [g-(y)FW)(y)]t}|0) , 



njf(,) f dxdye^''-''-''^'''-y 

.7 J 



X 



{Q\T{[q{x)T^'h'^\x)] Df\0) [q{y)r^'h(^\ym\0) . (3.16) 



In deriving Eq. (3.15) we have applied the relations 

I'J rfxd2/e-----'-^(0|T{[g(x)F^^/i(^)(x)]DfH0)[^^^^ (3.17) 

and 
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J dxdye-----'-^(0|T{[g(x)r^^/i(^)(x)]D;»(0 (3.18) 



Note that only the correlation function U.^^ is relevant to our purpose. It can be written in 
the double dispersion relation form 

The results of the QCD sum rules can be obtained in the following way. On the phe- 
nomenological side, which is the sum of the relevant hadron states, this correlation function 

can be written as 



where A is the binding energy of the heavy meson in the heavy quark limit and ellipses 
denote resonance contributions. On the theoretical side, the correlation function 11^'^ can be 
alternatively calculated in terms of quarks and gluons using the standard OPE technique. 
Then we equate the results on the phenomenological side with that on the theoretical side. 
However, since we are only interested in the properties of the ground state at hand, e.g., the 
B meson, we shall assume that contributions from excited states (on the phenomenological 
side) are approximated by the spectral density on the theoretical side of the sum rule, which 
starts from some thresholds (say, Uij in this study). To further improve the final result under 
consideration, we apply the Borel transform to both external variables lu and lu'. After the 
Borel transform [30] , 

BKU-.-')]= ^im Jim ;^(--')-"M^]-(--)'^^^[|;]'^nJf«(a;,.0, (3-21) 

■' —u>'—*oo — w— >oo ■ ■ ^ 

^ fixed ^ 

mt' 

the sum rule gives 

16 JO Jo 

where (Vij is the threshold of the excited states and pQ^^ is the spectral density on the 
theoretical side of the sum rule. Because the Borel windows are symmetric in variables ti 
and t2, it is natural to choose ti = t2- However, unlike the case of the normalization of the 
Isgur-Wise function at zero recoil, where the Borel mass is approximately twice as large as 
that in the corresponding two-point sum rule [31], in the present case of the three-point sum 
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rule at hand, we find that the working Borel windows can be chosen as the same as that in the 
two-point sum rule since in our analysis the output results depend weakly on the Borel mass. 
Therefore, we choose ^1=^2 = t. By the rcnormalization group technique, the logarithmic 
dependence asln{2t/ii) can be summed over to produce a factor like [as{iJ,) / as{2t)]'^ . After 
some manipulation we obtain the sum rule results: 



where 



16 



'asm' 



25 



-9/(2/3o) 



OPE 



i OPE,,,, 



1, i OPEs,, + OPE, 



(3.23) 



OPEf 



^(OPE) 



2 

2pt;i,j ) 



OPE 



^1, 



1 

16 



{qgsa ■ Gq) 



OPE 



0{as) , 



87r2 



iDTT"^ 



(3.24) 



with 



(OPE) 



2pt;i,j 



7r2 



ds s e 

2a 



as / 17 47r^ 
1 + — (— + 

TT V 3 9 



- 21n- 

t 



\ TT ' 



+ 



{qgsO_^Gq)_-) 
16^2 J ■ 



(3.25) 



For reason of consistency, in the following numerical analysis we will neglect the finite part 
of radiative one loop corrections in OPE^. ^ and OPE,. ^. (and in Eq. (3.28)). The parameter 
S in (3.23) is some combination of the P functions and anomalous dimensions (see Eq. (4.2) 
of [32]) and is numerically equal to —0.23. The relevant parameters normalized at the scale 
t are related to those at 11 by [33,31,35] 



F(2t)=F(/x)( 



asiiJi) ' 



a,(2t)s-2//3ol-5^ 



1 _ (5^£M ' 



a,(2t)N-4//3o 



OLs{[l) 

{gsQCT ■ Gq)2t = {gsqcr -Gq)^- ( — -— j 

{asG% = {asG^)^, 



(3.26) 
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where (• • •) stands for (0| • • • |0). In the calculation of the correlation function, we have 
also used the fixed-point gauge (the Fock-Schwinger gauge) x^A^{x) = with ^4^ being an 
external gluon field. Under this gauge, the generalized quark propagator in the external 
gluon field reads 



A 9sX -7 7'^ -7 

4 p-ruq p-rriq 



+ ..., (3.27) 

where a and h are the color indices, i and j the Lorcntz indices. 

Let us explain the results obtained in Eqs. (3.23) and (3.24). OPE^. is obtained by sub- 
stituting Dj^*^ by in n^f(i) [cf. Eq. (3.16)] and it can be approximately factorized as the 

3 

product of (0PE)2pt;ij with itself, which is the same as the theoretical part in the two-point 
F{ij) sum rule [31-33]. In the series of (0PE)2pt;i.j, we have neglected the contribution pro- 
portional to {qq)"^. (More precisely, it is equal to as(^q')^7r/324; see Ref. [31].) Nevertheless, 
the result of (OPE) b. in Eq. (3.24) is reliable up to dimension six, as the contributions from 
the {qqY terms in (0PE)2pt;ij are much smaller than the term (1 -|- 0:5/77)^(^5)^/16 that 
we have kept [see Eq. (3.25)]. Note that in (OPE)^. the contribution involving the gluon 
condensate is proportional to the light quark mass and hence can be neglected. Likewise, 
OPEg. is the theoretical side of the sum rule, and it is obtained by substituting Z^J*-*^ by 
T'" in Eq. (3.16). To the order of dimension-five, the main contributions to OPE^. are de- 
picted in Fig. 1. Here we have neglected the dimension-6 four-quark condensate of the type 
{qVX^q qV\"-q). It's contribution is much less than that from dimension-five or dimension- 
four condensates and hence unimportant (see [34] for similar discussions). It should be 
emphasized that nonfactorizable contributions to the parameters Bi arise mainly from the 
Qv _ jiv operator mixing. 

At this point, it is useful to compare our analysis with the similar QCD sum rule studies 
in [34] and [15]. First, Chernyak [34] used the chiral interpolating current for the B meson, 
so that all light quark fields in his correlators are purely left-handed. As a result, there 
are no quark-gluon mixed condensates as these require the presence of both left- and right- 
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handed light quark fields. The gluon condensate contribution enters into the £i sum rule 
with an additional factor of 4 in comparison with ours. Second, our results for OPE^. are 
very different from that obtained by Baek et al. [15]. The reason is that they calculated the 
full n^"" (obtained by replacing D^^^^ by in Eq. (3.14)) rather than the pseudoscalar part 
of n^*'". Therefore, their results are mixed with the 1+ to 1+ transitions. Also a subtraction 
of the contribution from excited states is not carried out in [15] for the three-point correlation 
function, though it is justified to do so for two-point correlation functions. Indeed, in the 
following analysis, one will find that after subtracting the contribution from excited states, 
the contributions of OPE^. are largely suppressed. Furthermore, as in the study of the B 
meson decay constant [31], we find that the renormalization-group effects are very important 
in the sum rule analysis. Consequently, there is no much difference between the resulting 
values of ei and £2- Moreover, Si aX, 11 — rrib are largely enhanced by renormalization-group 
effects. 

The value of F in Eq. (3.23) can be substituted by 



-A/t 



1 _ 25^2M 



26 



TT 



TT^ 



1 + ^ M h 21n-) 

. TT V 3 ^ 9 



/ 2a,(2t)x {qg,a ■ Gq)2t I „„s 
-(1 + —^) {qq)2t + 1 , (3.28) 



which can be obtained from the two-point sum rule approach [32,33,31]. For the numerical 
analysis, we use the following values of parameters [35,36] 

{qq),.=i GeV = -(240 MeV)^ 
{asG^)^.=l Gev = 0.0377 GeV^ , 

{q9s(Ti,.G^''q)^=i GeV = (0.8 GeV^) x cev , (3.29) 

as input and neglect the finite part of radiative one loop corrections. Since in our convention 
D/j, — df^ — igsA^, we have {gsqc ■ Gq) — ml{qq). Next, in order to determine the thresholds 
LOij we employ the B meson decay constant fs — (185 ± 25 ± 17) MeV obtained from a 
recent lattice-QCD calculation [37] and the relation [38] 

_2<m^., _(0.8~l.l)GeV. 

y/niB ^ 3 TT ^ ^ nib ' 

that takes into account QCD and l/m^ corrections. Using the relation between F{mb) and 
F{pl) given by Eq. (3.26) and = (4.85 ± 0.25) GeV, we obtain 
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F(ii = 1 GeV) ^ (0.34 ~ 0.52) GeV^/^ . 



(3.31) 



Since the A parameter in Eq. (3.28) can be replaced by the A sum rule obtained by applying 
the differential operator t^d\n/dt to both sides of Eq. (3.28), the -F(/x) sum rule can be 
rewritten as 

d 

F2(^) = (right hand side of Eq. (3.28)) x exp[t — ln(right hand side of Eq. (3.28))] , (3.32) 

which is A-free. Then using the result (3.31) as input, the threshold loq in the F{jji) sum 
rule, Eq. (3.32), is determined. The result for loq is 1.25 — 1.65 GeV. A larger = 1 GeV) 
corresponds to a larger uoq. The working Borel window lies in the region 0.6 GeV < t < 
1 GeV, which turns out to be a reasonable choice [33]. Substituting the value of luq back into 
the A sum rule, we obtain A = 0.48 - 0.76 GeV in the Borel window 0.6 GeV <t <1 GeV. 
This result is consistent with the choice rrib — (4.85 ± 0.25) GeV, recalling that in the heavy 
quark limit, A = rriB — Trih- To extract the cij''*'' sum rules, one can take the ratio of Eq. (3.28) 
and Eq. (3.23) to eliminate the contribution of F^/exp(A/t). This means one has chosen the 
same A both in Eq. (3.28) and Eq. (3.23). Since quark-hadron duality is the basic assumption 
in the QCD sum rule approach, we expect that the same result of A also can be obtained 
using the A sum rules derived from Eq. (3.23) (see [39,35] for a further discussion). This 
property can help us to determine consistently the threshold in 3-point sum rule, Eq. (3.23). 
Therefore, we can apply the differential operator t^dln /dt to both sides of Eq. (3.23), the 
^v{i) g^j^ Y^Ae^ to obtain new A sum rules. The requirement of producing a reasonable value 
for A, say 0.48 — 0.76 GeV, provides severe constraints on the choices of Ui^j. With a careful 
study, we find that the best choice in our analysis is 

a;. 1 = -0.02 GeV + ujq , a;i,2 = -0.5 GeV + ujq , a;2,2 = -0.22 GeV + ujq . (3.33) 

Applying the above relations with ujq — (1.25 ~ 1.65) GeV and substituting F(/x) in 
Eq. (3.23) by (3.28), we study numerically the (f^^^ sum rules. In Fig. 2, we plot and 

as a function t, where = 8df> /9 + 2df'^ /3, and e''i = -Adf^ /27 + 8df'^ /9. The dashed 
and solid curves stand for B^ and s^, respectively, where we have used ujq = 1.4 GeV (the 
corresponding decay constant is fe = 175 ~ 195 MeV or F(//) = 0.405 ± 0.005 GeV^/^). 
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The final results for the hadronic parameters Bi and are (see Fig. 2) ^ 

B'^ifx = 1 GeV) = 0.60 ± 0.02, B^{ii = 1 GeV) = 0.61 ± 0.01, 
£^(/x = 1 GeV) = -0.08 ± 0.01, £^(// = 1 GeV) = -0.024 ± 0.006. (3.34) 

The numerical errors come mainly from the uncertainty of lvq = 1.25 ~ 1.65 GeV. Some 
intrinsic errors of the sum rule approach, say quark-hadron duality or ag corrections, will 
not be considered here. 

Substituting the above results into Eq. (3.7) yields 

Bi{mb) = 0.96 ± 0.04 + 0{l/mb) , B2{mb) = 0.95 ± 0.02 + 0{l/mb) , 
ei{mb) ^ -0.U±0m + O{l/mb) , £2(^5) = -0.08 ± 0.01 + ©(l/mt) . (3.35) 



It follows from Eq. (2.11) that 
r{B-) 



1.11 ±0.02, 
1, 



= 0.99 



B 



) (0.007 + 0.020 S)r, 



(3.36) 



r{Bd) ^185 MeV' 

to the order of l/m\. Note that we have neglected the corrections of SU(3) symmetry break- 
ing to the nonspectator effects in t{Bs)/t{B(i). We see that the prediction for t{B~) /t^B^) 
is in agreement with the current world average: T{B~)/T{Ba) =1.07± 0.04 [7], whereas 
the heavy-quark-expansion-based result for t{Bs)/t{B(i) deviates somewhat from the cen- 
tral value of the world average ^: 0.95 ± 0.05. Thus it is urgent to carry out more precise 
measurements of the Bg hfetime. Using the existing sum rule estimate for the parameter 
r [10] together with B = 1 gives T{Ai))/T{Bd) > 0.98. Therefore, the 1/ml nonspectator 
corrections are not responsible for the observed lifetime difference between the and B^. 



^For comparison, the sum rule results obtained in [15] and [34] are eY^^^{fi) = —0.041 it 0.022, 
efLLS^^) = 0.061 ± 0.035 [15] and (/Lt) « -0.15, £^{n) « [34], respectively, with being a 
typical hadronic scale ~ 0.70 GeV. Note that the definition of -Bi(/i) and £i(/x) in [15,34] is different 
from ours by a factor of F"^ {rrib) / F"^ (fj,) , that is, ef^^^'^{n) = £,(//) x F'^ {rrib) / F'^ {n) . 

^For example, the neutral B meson lifetimes are measured at CDF to be [40]: T{Bd) = 1.58 =b 
0.09 ± 0.02 ps and t{Bs) = 1.34l°;23 ± o.05 ps. 
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IV. DISCUSSIONS AND CONCLUSIONS 



The prediction of B meson lifetime ratios depends on the nonspectator effects of order 
Idn^/ml in the heavy quark expansion. These effects can be parametrized in terms of the 
hadronic parameters Bi, B2, £1 and 62, where Bi and B2 characterize the matrix elements of 
color singlet-singlet four-quark operators and £1 and £2 the matrix elements of color octet- 
octet operators. 

As emphasized in [12], one should not be contented with the agreement between theory 
and experiment for the hfetime ratio t{B^)/t{B(i). In order to test the OPE approach for 
inclusive nonleptonic decay, it is even more important to calculate the absolute decay widths 
of the B mesons and compare them with the data. From (2.6) and (3.35) and considering 
the contributions of the nonspectator effects, we obtain 



noting that the next-to-leading QCD radiative correction to the inclusive decay width has 
been included. The absolute decay widths strongly depend on the value of the h quark mass. 

The problem with the absolute decay width V{B) is intimately related to the B meson 
semileptonic branching ratio i3sL- Unlike the semileptonic decays, the heavy quark expansion 
in inclusive nonleptonic decay is a priori not justified due to the absence of an analytic 
continuation into the complex plane and hence local duality has to be invoked in order to 
apply the OPE directly in the physical region. If the shorter hfetime of the A;, relative to 
that of the B^, meson is confirmed in the future and/or if the lifetime ratio t{Bs) /T{Bd) is 
observed to be approximately unity, then it is very likely that local quark-hadron duality is 
violated in nonleptonic decays. It should be stressed that local duality is exact in the heavy 
quark limit, but its systematic l/mq expansion is still lacking. Empirically, it has been 
suggested in [9] that the presence of hnear l/rri}, correction, described by the ansatz that the 
h quark mass m;, is replaced by the decaying bottom hadron mass in the factor in 
front of all nonleptonic widths, will account for the observed lifetime difference between the 
Aft and B^- To be specific, the ansatz Fnl ^nhijUHi^/mb)^ will lead to the results [12]: 



Ttoi{Ba)= (3.6llJ:°l)x 10-^3 GeV, 
rtot(S-)= (3.34lJ:°l)xlO-i^GeV, 



(4.1) 



^(A.) 



0.76, 




0.94. 



(4.2) 
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This simple prescription not only solves the lifetime ratio problem but also provides the 
correct absolute decay widths for the and the B mesons. The predicted lifetime hierarchy 



r(A,) > t{E;) > r(SO) > r(Q,) 



(4.3) 



for bottom baryons is in sharp contrast to the OPE-based hfetime pattern [12]: 



r(Q,):^r(S,-)>r(A,)^r(S°). 



(4.4) 



Of course, whether this empirical ansatz truly works or whether it can be justified in a more 
fundamental way (see, for example, [42]) remains to be investigated. Nevertheless, it is worth 
emphasizing that, although a linear l/mg correction to the inclusive nonleptonic decay rate 
is possible [43,44], the violation of local quark-hadron duality docs not necessarily imply the 
presence of I/ttlq terms in inclusive widths and hence the aforementioned ansatz. 

To conclude, we have derived in heavy quark effective theory the renormalization-group 
improved sum rules for the hadronic parameters Bi, B2, £1, and £2 appearing in the matrix 
element of four-quark operators. The results are Bi{mb) = 0.96±0.04, B2{mb) = 0.95±0.02, 
£i{nib) = — 0.14±0.01 and £2{nib) = — 0.08±0.01 to the zeroth order in 1/mf,. The resultant 
5-meson lifetime ratios are T{B~)/T{Bd) = 1.11 ± 0.02 and T{Bs)/T{Bd) ~ 1. 
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(a) (b) 

FIG. 1. The main diagrams contributing to OPEe- [cf. Eq. (3.24)]: (a) the contribution from 
the gluon condensate, and (b) the contribution from the quark-gluon mixed condensate. In (b) the 
mirror-symmetric diagram is included in the calculation. The double lines denote heavy quarks in 
HQET. 
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FIG. 2. and £V(;x) as a function t, where BV = + 2df^ and 

= + 8cZ2^'V9- The dashed and soUd curves stand for B'V and e", respectively. Here 

have used ojq = 1.4 GeV and Eq. (3.33). 



24 



